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Abstract. In this paper we report an experimental investigation of various statis- 
tical properties of the spatial Fourier modes of the vorticity field in turbulent jets 
for a large range of Reynolds numbers (530 < R\ < 6100). The continuous time 
evolution of a spatial Fourier mode of the vorticity distribution, characterized by 
a well defined wavevector, is obtained from acoustic scattering measurements. The 
spatial enstrophy spectrum, as a function of the spatial wave- vector, is determined by 
scanning the incoming sound frequencies. Time- frequency analysis of the turbulent 
vorticity fluctuations is also performed for different length scales of the flows. Vortic- 
ity time-correlations show that the characteristic time of a Fourier mode behaves as 
the sweeping time. Finally, we report preliminary Lagrangian velocity measurements 
obtained using acoustic scattering by soap bubbles inflated with helium. Gathering 
a very large number of passages of isolated bubbles in the scattering volume, one 
is able to compute the Lagrangian velocity PDF and velocity spectrum. Despite 
the spatial filtering due to the finite size of the bubble, the latter exhibits a power 
law with the -2 exponent predicted by the Komogorov theory, over one decade of 
frequencies. 
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1. Introduction 

In fully developped turbulence, most of the experimental, numerical 
and theoretical works [1] rely on the statistics of the longitudinal ve- 
locity increments to study the dynamic of a given turbulent length 
scales in the physical space. In the Fourier space, there are numerous 
theoretical and some numerical studies ([2]), but very few experiments. 
The objective of this paper is mainly to give experimental vorticity data 
in the Fourier space. Firstly, the ultra-sound acoustic scattering method 
is described, and the experimental conditions are detailed. Then, both 
time and spatial spectra of vorticity are presented and are extended to a 
time-frequency analysis. Correlation time of the vorticity Fourier mode 
is discussed. Finally, we show how the same scattering method can be 
applied to detect and follow the trajectory of isolated soap bubbles in 
a turbulent air jet, in order to extract Lagrangian velocity data and 
compute statistics. 

© 2008 Kluwer Academic Publishers. Printed in the Netherlands. 
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1.1. Acoustic scattering by vorticity 

Wave propagation (light and sound) in fluids is known to be strongly 
affected by turbulent velocity and/or temperature gradients. In the 
case of direct propagation, turbulence is responsible for random wave 
amplitude and phase fluctuations along the ray paths leading to the well 
known stellar scintillation phenomena limiting the spatial resolution 
of optical and radio observations [2, 3]. Acoustic wave propagation is 
also sensitive to local velocity and temperature fluctuations induced by 
turbulence: such fluctuations are responsible for local fluctuations of 
the sound velocity along the ray paths which result in distorsions of 
the wave fronts [2, 4, 5, 6]. 

In the presence of spatial inhomegeneities with length scales compara- 
ble to the acoustic wave-length \ s , an incident acoustic plane wave can 
also be scattered giving rise to scattered acoustic waves propagating in 
directions away from the incident wave direction of propagation [7, 8]. 
There exists a large number of theoretical [9, 10, 11, 12, 13] and numeri- 
cal [14, 15, 16] studies dealing with the acoustic scattering phenomenon 
by velocity fields. Most of it was initiated nearly fifty years ago and 
started with the papers authored by A.M. Obukhov [9], R.H. Kraich- 
nan [10], B. T.Chu and L. S. G. Kovasznay [11], G. K. Batchelor [12] 
to mention but a few. Recently, F. Lund [13] has established, under 
reasonable and fairly restrictive assumptions, a linear relation between 
the scattered amplitude of a plane acoustic wave incident on a turbulent 
flow and the spatial Fourier transform of the vorticity field. 
In analogy with the more usual light scattering phenomenon, and fol- 
lowing Batchelor [12] the physical mechanism at the origin of acoustic 
scattering by vorticity can be thought of as follows: an acoustic wave 
impinging on a vorticity distribution induces fluctuations of the vortic- 
ity at the incoming sound frequency (by virtue of the Kelvin circulation 
theorem) . Each scatterer (vortex) acting as a secondary source will, in 
turn, radiate a sound wave. The coherent average (taking into account 
the relativepositions of the individual vorticity elements) over the scat- 
terers distribution results in the emission, outside the vorticity domain, 
of a scattered acoustic waves. Note that whereas the light scattering 
process is usually linear, the acoustic scattering phenomenon depicted 
here stems from the non-linear term of the Navier-Stokes equation, 
and requires more detailed explanations and computations, beyond the 
scope of this paper. In particular, one can find in [13] a clear explanation 
of the respective contributions to the total acoustic scattering ampli- 
tude of the vorticity field on one hand, and of the irrotational velocity 
field, induced by the vorticity field, on the other. Using a Born approx- 
imation, Lund et al. obtain the following linear relation between the 
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scattered acoustic pressure amplitude and the spatial Fourier transform 
of the vorticity field: 



Pscat(v) 2 . -COs(9 scat ) y^D/c 

Pscat^ =7r 2 Kscat)_ — (n A r).fl(q sca t, V — v ) (1) 

Pine l-COS[B scat ) C l U 

where A and "." stand for the vector product and the scalar product 
respectively, and where the scattering wave- vector (also called momen- 
tum transfert) is given by: 



27T Attv . 9 scat r - n 

Iscat = — {vr - u n) ~ S m(— — )- for v ~ v Q (2) 

c c 2 r — n 




Figure 1. Acoustic Scattering by a Vorticity Distribution 

In the above equations, the meaning of the different variables is the 
following according to the schematic drawing in fig. (1): 

— <iscat = k — k D is the scattering wave-vector. 

— n and r are the unit vectors in the incident and scattered wave 
directions of propagation, respectively. 

— k Q (resp. v a ) is the vector wave-number (resp. frequency) of the 
incoming sound wave (in the direction n). 

— k (resp. v) is the vector wave-number (resp. the frequency) of the 
scattered sound wave (in the direction of observation r) . 

— sca t is the scattering angle. 
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— pi nc is the complex pressure amplitude of the incoming sound wave 
(assumed to be plane and monochromatic). 

— Pscat is the complex pressure amplitude of the scattered sound 
wave. 

— D, c and A s stand respectively for the acoustical path between the 
measurement area and the detector, the adiabatic sound velocity 
and the acoustic wavelength. 

1.2. Acoustic scattering: a probe of vorticity 

Equation (1), indicates that, using an incident acoustic plane wave 
(ideally of infinite extension), one can directly probe the time-space 
Fourier mode f2j_(q sca t, / = v — u Q ) of a well defined component of the 
vorticity field fij_(r, t): 

Ol(<W,/) = /// Jn ± (r,t)e-^* r - 2 ^dtd 3 r (3) 

Actually, the direction of the probed vorticity component (indicated 
by the subscript j_) is perpendicular to the scattering plane defined 
by the vector wave-numbers of the incident and scattered (detected) 
acoustic waves. Equation (1) can be reformulated as a relation between 
the scattered amplitude and the convolution product of the incident 
acoustic amplitude and fi_L(q SC atj f = v — u Q ). Using an inverse Fourier 
transform on the frequency variable, one then obtains the following 
linear relation between the time variables : 

p.scat(K t) oc L(e scat )n ± (q scat = k - k D , t).p inc (k , t) (4) 

expressing the scattered pressure amplitude as the result of a mod- 
ulation of the incident acoustic pressure by the space Fourier trans- 
form f2(q sca t,i). The variations of the angular prefactor L(9 scat ) with 
the scattering angle 9 sca t, represented on figure (2), are typical of a 
quadrupolar like radiation pattern, diverging at null scattering an- 
gle (where the Born approximation breaks down), and exhibiting null 
values at 9 scat = 90° and 9 scat = 180° (back-scattering). 

Among the several hypothesis put forward in [11, 13] to arrive at 
Eqn. (4), necessary for the Born approximation to be applicable, the 
velocity amplitude u soun d = of the incident wave and the Mach 

number M = ^H^- of the flow under investigation (associated to the 
probed vorticity distribution) must be such that : u soun d <C Vfi ow <C c. 
In addition, the time scales of the flow are required to be much smaller 
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than the period T Q = of the incident wave. These conditions can be 
easily fulfilled using ultrasonic acoustic waves with low enough inten- 
sity. From a practical point of view, equations (1) and (4) imply that 
the scattered pressure signal is a narrow-band signal centered around 
the incident frequency v Q \ the continuous time evolution of f2_i_(q sca t, t) 
can then be easily recovered by a simple demodulation operation using 
for exemple an heterodyne detection (either analog or digital) [17]. 
After demodulation, for a fixed incoming sound frequency v Q and a 
fixed scattering angle 6 SC at one thus gets the time dynamics of a spatial 
Fourier mode of the vorticity tt± (q SC(l t , t) , caracterized by a well defined 
vector wave- number q scat (with specified modulus and orientation). 
The overall effect of the acoustic scattering process can then be viewed 
as a time-continuous spatial band-pass filtering operation of the vortic- 
ity distribution Q±(r,t). The center frequency of the spatial filtering 
operation is given by the scattering vector wave-number q sca t which 
depends on a combination of the incident frequency v Q and of the 
scattering angle 9 SC at as stated by Eqn (2). For practical reasons and 
as we dispose of wide-band transducers (emission and reception), the 
analysed length-scale of the flow is tuned through the selection of an 
incident frequency v Q while the scattering angle is usually kept con- 
stant. It is worth noting here that, as a result of the demodulation 
process according to which the instantaneous scattered pressure signal 
is multiplied by both the in phase (cos{2i:v t)) and the in quadrature 
(sin(2Trv t)) signals phase-locked with the electric signal driving the 
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sound transmitter we get finally a low frequency signal which is com- 
plex (modulus and phase) and directly proportionnal to the complex 
quantity £l±(q scat , t), where: 



Thus, the instantaneous phase of the demodulated scattered pres- 
sure signal is simply a measure of the phase modulation of the scattered 
pressure wave with respect to the incident wave (chosen as a phase 
reference). Similarly, the instantaneous amplitude of the demodulated 
scattered pressure signal is a measure of the amplitude modulation 
with respect to the constant amplitude of the incident wave. Both in- 
formations (phase and amplitude) are easily recovered from the electric 
signals owing to the linearity property of the acoustic transducers used 
to generate and detect acoustic waves. 

1.3. Spatial filtering: finite size effects and spectral 
resolution 

As indicated by equation (1) and figure (2), the acoustic amplitude 
scattered by vorticity distributions is a strongly decreasing function of 
the scattering angle 9 scat . In order to preserve an acceptable signal-to- 
noise ratio, one needs to work with low scattering angles (typ. between 
10° and 60°). Hence, we work with a bistatic configuration, using a 
pair of acoustic transducers : one transmitter and one receiver. The 
relative positions of the two tranducers, as well as their arrangement 
with respect to the flow under investigation define the scattering angle 
and the orientation of the scattering wave-vector c[ sca t with respect 
to the flow geometry (e.g. its mean velocity). A typical experimental 
configuration, for the study of a turbulent jet flow is represented on 
figure(3). Equation (1) is valid for an ideal incident plane wave and a 
perfectly directive detection. In a real experiment, transducers are of 
finite dimensions and one has to take into accountdiffraction effects. Let 
L be the typical size of the acoustic transducers and X s = — the acous- 
tic wavelength, at large distances from the transducer (in the far-field 
limit D > ^ [18, 4]), one expects acoustic energy to be either emitted 
by the transmitter or detected by the receiver in a cone corresponding 
to a continuous distribution of angles defining the angular aperture 
(or directivity). In particular, the principal diffraction lobe, around the 
normal to each tranducer, is a cone with a half angle aperture A8 
given by [18, 4]): A6 ~ sin(A6) = By differentiating Eqn. (2), one 
can associate to these distributions of emission and detection angles, a 
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distribution of probed scattering wave- vectors q sca t = 4n — with 

a typical width: 

Aq scat = 2vr^A0 ~ (6) 

As — TW 
cos(-f 

Equation (6) gives an estimation of the width of spatial band-pass 
filtering operation performed by the scattering process around the spa- 
tial frequency (length scale) q sca t. It is interesting at this point to 
notice that Equation (6) can be interpreted as a mere consequence of 
the uncertainty relation Aq sca tAx ~ 2ir, where Ax = g (fj ^ is roughly 



the extension of the measurement volume (in the direction of q SC at)> 
defined by the intersection of the incident and detected (antenna beam) 
acoustic beams (cf fig. (3)). An equivalent formulation, albeit in the 
physical space, of the previous result can be found for example in Ref. 
([4] pp. 443-446), where the scattering amplitude in an experimental 
bistatic configuration is expressed as: 



r r r+oo 

Ol(<W,*) = J J J \F rec (r)F trans (r)\n ± (r,t)e-^ r d 3 r (7) 

where G(r) = \F rec (r)(F trans (r)\ is a ponderation function defined 
as the modulus of the cross product of the complex amplitudes of the 
incident (resp. detected) F rec (r) (resp. F trans (r)) acoustic beams. 

Let's note V sca t, the volume of the region of space over which G(r) 
takes appreciable non zero values and let's call V sca t the scattering 
volume, equation (7) can be expressed as a windowed spatial Fourier 
transform: 

Ol(<W,*)= / / / n ± (r,t)e-^ T d 3 r (8) 

J J J Vscat 

In conclusion, we deal with a complex signal z scat (t) = p{t)eP^^ oc 
Q±(ciscat,t) with a spectral resolution given by : 

Qscat = Qscat^x ± Aq x e x ± Aq y e y ± Aq z e z (9) 
qscat = —^-sin(9 scat /2) (10) 

Aq x ~ Aq y ~ Aq z ~ (t^—) 1/3 ^ \ (H) 

Vscat cos(e scat /2) 

In the following, we shall refer to p sca t(t) as the scattered pres- 
sure signal (real) with frequencies around v (high frequency) and to 
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Figure 3. Acoustic scattering experiment in a turbulent jet flow 



Zscat(t) as the demodulated scattered pressure signal (complex) with 
low frequencies. 

2. Experimental Conditions 

2.1. Flow configurations 

Two axisymetric turbulent jets have been investigated: a laboratory 
air jet at moderate Reynolds number, (hereafter, called LEGI jet), and 
a low temperature gaseous helium jet at very high Reynolds number 
which has been performed in the cryogenic facility hosted by the CERN 
in Geneva. In the laboratory air jet the center of the measurement 
volume was located on the axis at about 40 diameters downstream the 
nozzle (with a diameter of 12 cm). The mean velocity and the rms veloc- 
ity have been measured in the center of scattering measurement volume, 
with conventional instrumentation (hot wire and Pitot anemometry). 
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Table I. Turbulent flows characteristics 



Flow 


Rx 


A (mm) 




to 
L 


27T J Cfmin 


2"7T 1 Qmax 


LEGI Jet 


530 


8.4 


180 


0.75 


~ 3.5A 


~ 10r7 




740 


6.7 


120 


11 


~ 4.5A 


~ I5r) 




785 


6.5 


107 


11 


~ 5A 


~ 17rj 


CERN Jet 


3450 


0.69 


5.45 


1,5 


~ 10A 


~ 1.2A 




4750 


0.51 


3.40 


11 


~ 15A 


~ 1.8A 




6090 


0.42 


2.55 


11 


~ 17A 


~ 2A 



The value of the Taylor microscale Reynolds number estimated from 
longitudinal velocity samples and isotropic and homogeneous relations, 
ranges 530 < R x < 785 [19]. 

The high Reynolds jet investigated with acoustic scattering [20, 21] 
consists in an axisymetric gaseous helium jet at a temperature of about 
A°K. The flow emerges from a 2.5 cm diameter nozzle into a large 
thermally insulated cylindrical chamber (4.6 m high and with a diam- 
eter of 1.4 m). The cryostat is connected to a large refrigerator, with 
6 kW cooling power, allowing flow rates of cryogenic Helium as high 
as 250 g/s. Taking advantage of the very low kinematical viscosity of 
Helium at such a low temperature (v ~ 810 -8 m? / s at T = 4.5 K 
and P = 1.2 bar), very large Reynolds numbers can be achieved (up to 
R\ ~ 6000) . Longitudinal velocity fluctuations have also been measured 
with a dedicated super-conducting 'hot-wire' anemometer located at 
the center of the acoustic measurement volume, on the jet axis, 50 
nozzle diameters downstream. For the highest flow rate the mean ve- 
locity in the measurement volume was ~ 4.4 m/s and the rms velocity 
~ 1.25 m/s. 

The main flow characteristics of the different experiments are summa- 
rized in table I, as well as the range of investigated spatial wave-vectors. 

2.2. Experimental Set-Up 

Acoustic waves are generated and detected using Sell-type transducers 
(see e.g. [22]). The Sell tranducers are electro-acoustic reciprocal trans- 
ducers, consisting in a circular plane piston, of diameter L = 15 cm 
and L = 5 cm, respectively in the LEGI and the CERN jet flows, 
polarized with a high static voltage (« 200 Volts DC). The circu- 
lar membrane made of a very thin (thickness ~ 15 fim) mylar sheet 
achieves a large frequency band-width (between 1 kHz and 200 kHz in 
air flows). Indeed, when the working fluid is air, frequencies higher than 
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200kHz are strongly attenuated by dissipative process (viscosity and 
thermal conductivity) . At small acoustic frequencies (below 1 kHz) , the 
main frequency limitations stem from the turbulent noise generated by 
the turbulent flow under investigation and the beams divergency due 
to diffraction effects. The transmitter is driven by a MOSFET power 
amplifier (NF 4005, 100 VA), while the acoustic pressure signal collected 
on the receiver is converted into a voltage signal using a home made 
linear, low noise charge amplifier. 

In all experiments, the scattering angle 9 sca t is kept at a constant value 
(typ. 60° in air flows) chosen so as to realize a trade-off between op- 
timizing the sensitivity of the measurements and limiting the spurious 
effects induced by diffraction side-lobes. Moreover, in any case, the size 
L of the transducers is such that the linear size (T^cat) 1 ^ 3 is of the 
same order than the integral length scale of the turbulent flow (note 
that in the CERN experiment, the sound velocity is about 110 m/s, 
nearly three times smaller than in the air at usual temperature). In a 
typical sound scattering experiment time series of the acoustic pressure 
signal are collected by the Sell receiver and then sampled and recorded 
using Agilent E1430 VXIbus-based analog to digital converters with 
a high precision (23 bits). Each E1430 module is provided with an 
analog anti-aliasing filter, digital filtering and decimation circuits and 
its own local oscillator allowing real-time heterodyne demodulation. In 
order to reduce phase noise and frequency drift, the sampling clock, 
the local oscillator clocks of the E1430 digitizer and the clock of the 
frequency generator (Agilent E33120) which generates the continuous 
sinusoid signal driving the transmitter are locked on the same 10 MHz 
master clock. The frequency of the local oscillator used for the numer- 
ical demodulation by the E1430 digitizer is thus precisely tuned to the 
frequency v Q of the incoming sound wave. 



3. Spectral Analysis and Spatial Enstrophy Spectrum 

3.1. Temporal and Spectral Caracteristics of the 
Scattered Pressure Signal 

We consider in this section, the averaged statistical spectral proper- 
ties of the demodulated scattered pressure signal, recorded on the 
detector for fixed incoming sound frequency v Q and scattering angle 
9 sca t, defining a single analysing wave- vector <\ scat according to Eqns. 
(1, 2). Indeed, being a direct image of the time evolution of Q±(r,t) 
which is a stochastic signal, one expects the demodulated scattered 
pressure signal z scat (t) to be also a complex random signal. We can 
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express the complex signal z scat (t) as p(t)e l< ^*), where p(t) (resp. (f>(t)) 
is the instantaneous amplitude modulation (resp. phase modulation). 
Typical evolutions of these two quantities are represented on figure (4) , 
obtained in the LEGI jet for a scattering wave- vector (q sca t^ = 4.93). 
The upper plot shows the evolution, along time, of the instantaneous in- 
tensity I(t) = \z sca t\ 2 = p{t) 2 , proportional to the instantaneous spatial 
enstrophy |fi(q sca t, t)\ 2 which clearly exhibits very large fluctuations. 
The lower part of figure (4) sketches the time evolution of the instan- 
taneous phase 4>(t). In this representation, the phase 4>(t) has been 
unwrapped, using an appropriate algorithm, to take into account phase 
jumps greater than it. The behaviour of the phase shift is roughly linear 
with a mean slope = 5800 rad/s. Indeed, since in this experiment 
the direction of the scattering wave-vector q sca t has been chosen to be 
aligned with the mean flow velocity vector V ai , 9 , the linear phase shift 

is a mere consequence of a Doppler effect 5u{t) = v(t) — v Q = ^^jp- 
with : 

d( t>(t) _ v (m 

^ — Hscat v avg 

This Doppler effect is related to the advection of the vorticity dis- 
tribution by the large scale velocity field. Note however, that local 
fluctuations of the instantaneous frequency shift are still visible, owing 
to turbulent fluctuations in time of the local velocity itself. From the 
knowledge of q sca t (739 m~ l ) and the estimated slope of , we get 
an estimation of the local mean flow velocity : 7.79 m/s of the order 
of the mean flow velocity obtained with usual hot-wire measurements, 
averaged on a much longer time (7.39 m/s). 

Let's turn now to the spectral properties of the demodulated scat- 
tered pressure signal. From a long time serie (2 20 points sampled at 
F s = 16384 Hz, corresponding to 2500 integral times Ti nt ), we estimate 
the average Power Spectral Density (the spectrumhereafter will be re- 
ferred as PSD sca t{5v)) using the usual Welch's averaged periodogram 
method [23] (with a hanning window of length 2048 samples and an 
overlap between consecutive segments of 1024 samples). The resulting 
PSD estimation, for the previous signal is plotted on figure (5). The 
spectrum PSD sca t(5v), is asymmetric (with respect to the incident 
frequency v Q (corresponding to a null Doppler frequency bv = 0) and 
has a Gaussian shape centered on a positive Doppler shift frequency 
as we choose <\ sca tV avg > 0. Note that, thanks to complex nature of 
the demodulated signal z sca t(t), we are sensitive to the direction of the 
flow: we are thus able to discriminate positive (i.e. downstream) and 
negative (i.e. upstream) velocities. Indeed, the Gaussian shape of the 
scattered pressure PSD scat (5v) is a direct consequence of the Gaussian 
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Figure 4- Instantaneous scattered pressure intensity (upper) and phase modulation 
(lower) in the LEGI jet flow (R x ~ 785) 



statistics of the large scale velocity. Using a non linear Gaussian fit (the 
dashed line on figure (5)) : 



PSD scat {5v) 



A{v ) 



2Tr5u r 



-exp 



2(Sis r ms) 2 



(13) 



we find a value 5v. 



avg 



865.5 Hz and a value bv r 



224.5 Hz 



measuring the width of the fluctuations of the advection velocity around 
its mean value. The estimated values 5v avg and 8v rms can be con- 
verted into values of the advection velocity according to Eqn. (12) : 

Vr, 



avg 



7.37 rajs and V„ 



1.91 m/s in very good agreement 



with hot-wire anemometry measurements of the longitudinal velocity 
(7.39 m/s and 1.89 m/s). 
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Figure 5. Power Spectral Density of f2x(q SC at,£) at a scale ~ A in the LEGI jet 
(7? A ~ 785) 



3.2. Power spectral density of the vorticity Fourier 
modes in high reynold numbers turbulence 

Acoustic scattering measurements have been also conducted in the 
CERN jet. Indeed, taking advantage of the very low kinematical vis- 
cosity of gazeous Helium at low temperature (around 4 °K), it is 
possible to achieve large Reynolds numbers, in a well controled flow 
(implying statistical stationnarity of the turbulent flows under inves- 
tigation), while keeping moderate Mach number values in order for 
the Born approximation used in acoustic scattering analysis to apply 
Figure (6) displays the PSD of the scattered acoustic pressure signal for 
R\ = 3450). For this figure, the scattering angle is set to 9 scat = 30° 
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and the incoming sound frequency to v Q = 110 k Hz giving a scattering 
wave-number q sca t = 3250 mT 1 equivalent to an anylsing scale 2.8A. 



-50 




-8 -4 4 8 

Doppler Frequency 8v (kHz) 

Figure 6. Power Spectral Density of f2j_(q a cat, i), at a scale ~ 3 A in the CERN 
helium jet flow ^ 3450). 

As in the LEGI jet, a large amplitude scattered peak is still visible 
on the PSD; when fitted with a Gaussian, one gets an estimation of the 
mean advection velocity of 1.24 m/s, in agreement with longitudinal 
velocity statistics obtained by hot-wire anemometry. Notice, however, 
that from the variance of the Doppler shift frequency 5v rms we obtain 
an estimation of the variance of the velocity significantly larger than 
the one given by hot-wire anemometry. The observed discrepancy, of 
order 20%, could be ascribed to the finite size of the transducers. 
A salient feature of the PSD at large Reynolds number, lies in the 
presence of significative scattered pressure amplitude at high Doppler 
frequency shifts, manifesting as exponential tails on both side of the 
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principal Gaussian peak. Note that a close look at figure (5), also reveals 
a similar behavior, albeit with a much smaller amplitude (relative to 
the Gaussian peak). We have checked that the two exponential- like 
wings are approximately symetric with respect to the mean Doppler 
shift frequency, indicating that this high frequency part of the PSD is 
associated with vorticity events, advected by the large scale flow. As 
the typical frequencies involved in these exponential wings are an order 
of magnitude larger than the Doppler shifts, we are inclined to ascribe 
them to a much richer dynamic of the time evolution of the small 
scale vorticity events. Comparisons with PSD obtained at still higher 
Reynolds numbers in the CERN experiment (up to R\ = 6090), 
indicate a clear enhancement of this high frequency behaviour (with 
respect to the low frequency Gaussian part). 

3.3. Spatial Enstrophy Spectrum 

In the LEGI Jet experiment, the scattering measurements have been 
repeated for various incoming sound frequencies (the scattering angle 
being held to the constant value 60 deg) in order to scan a large range of 
spatial wave- vector numbers. For each incoming sound frequency v , the 
time spectrum PSD sca t(Sv) has been fitted with a Gaussian function 
(according to Eqn. (13)) in order to estimate Sv avg (q sca t), 5v rms (q scat ) 
and A(q scat ). Figure (7) shows a plot of the evolution of Sv aV g(q sca t) and 
8v rms (q sca t) with the wave- number q sca t- A linear evolution is roughly 
evidenced, according to Eqn. (12). 

By integrating the scattered pressure PSD with respect to the Doppler 
shift frequency bv, one gets for each wave- number q sca t, a measure of the 
spatial enstrophy spectrum according to the Wiener-Kintchine theorem 
[17] : 

+T/2 +oo 
T^Sof / \ n (lscat,t)\ 2 dt= J \tt{q scaU 5v)\ 2 d{5v) (14) 
-T/2 -oo 

Using Eqn. (1) and Eqn. (13): 

+oo 

J \z scat (6v)\ 2 d(5v) = A(v ) (15) 

— oo 

one sees that the quantity q SC at 2 A(u ) , is a direct estimation of 

the spatial enstrophy spectrum Enstro(q sca t) = qscat 2 (\Q(q S cat, t )\ 2 ) t 

at wave-number q scat = ^psh^^p) up to a multiplicative factor 
\H(u )\ 2 which accounts for the frequency response of both transmitter 
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Figure 7. Evolution of the mean 8v avg and of the variance 5v rma of the Doppler 
shifts, with the scattering wave-number q aca t in the LEGI jet (R\ ~ 785) 



and receiver around the working frequency v Q . Actually, the transfer 
function H(v ) is evaluated from measurements of the the electrical 
response delivered by the receiver when the transmitter is driven by 
a wide-band electrical white noise. We have plot on figure (8) the 
variations of Enstro(q sca t) against the turbulent wave- numbers q SC at, 
for various length scales in the inertial range (R\ = 785). 

This log-log representation suggests a power-law scaling, with an 
exponent slightly positive and close to Kolmogorov 1941 prediction 
(dashed line with a slope 1/3). Resorting to a scale invariance argu- 
ment and assuming the statistical homogenity and isotropy of the flow, 
one expects (following Kolmogorov) the spatial enstrophy spectrum to 
scale as : Enstro(k) tx k 2 E(k) where E(k) oc A;~ 5 / 3 is the turbulent 
energy spectral density. This latter result is the mere consequence of 
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Figure 8. Spatial enstrophy spectral density as a function of the non dimensionnal 
wave-number in the LEGI jet (R\ ~ 785). 



the definition of the vorticity field fi(r, t) = V A u(r, t) which leads to 
fi(k, t) = — A u(k, t) in the Fourier space. 

The spatial enstrophy spectrum differs significantly from time-spectra 
measured with multi- wires velocity probes [24, 25, 26]. In its principle, 
the acoustic scattering measurement gives a good resolution in the 
Fourier space, thus it gives the true ID spectrum, as the spatial Fourier 
transform of the vorticity field is evaluated at a well defined spatial 
wave-vector. In the case of multi-wires anemometry measurements of 
the Eulerian vorticity, the vorticity spectrum evaluated at wave-number 
q x = 2-Ktt— involves the contributions of all the spatial wave- vectors 

vavg 

with the ID projection q x [27, 28]. Taken as a whole, such a 'ID spec- 
trum' is mainly dominated by the small scales (large q x ) viscous cut-off 
and hence does not exhibit any inertial scaling [29]. 
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3.4. Time-Frequency Distributions 

A detailed examination of the time evolution of the scattered pressure 
signal (phase and amplitude as depicted on figure (4)), reveals strong 
fluctuations which can be ascribed to the passage in the measurement 
volume of finite duration vorticity events. In turn the scattered acoustic 
signal exhibits a time-varying spectrum, as the frequency of the mea- 
sured pressure signal drops from v$ to i/q + bv when a vorticity event is 
present, with significative energy at the wave- vector q sca t in the spatial 
Fourier domain. Due to the randomness of the advection velocity of 
the vorticity field, the true vorticity spectrum is blurred and wiped 
out by the time averaging process. Thus we seek after a representation 
of the signal which preserves simultaneously the time and frequency 
information. 

Time- Frequency distributions [30, 31, 32], among which the Wigner- 
Ville distribution (WVD) is the most famous, provide a simultaneous 
description (in time and frequency) of the energy distributions of non- 
stationary signals in the time-frequency plane. The Wigner-Ville dis- 
tribution W z (t, v) of a complex analytic signal z(t) is defined as the 
time Fourier transform of the local auto-correlation function of z(t): 
R z (t,r) = z(t + ^)z*(t — ^) with respect to the lag variable r: 

W z (t, v) = F T [z(t + T -)z*{t - T -)\ = F T [R z (t, r)\ (16) 

where z(t) is the time signal and z*(t) its complex conjugate. The 
operator Fr\ denotes the Fourier transform operator, with respect to 
the variable (.). This definition generalizes, in some way, the Wiener- 
Kinchine theorem to non-stationary signals, for which the auto-correlation 
function depends on two time variables: R z (h,t2) = {z(t\)z* (£2)) / 
Rz(ti—t2)- From an experimental and heuristic point of view, the effect 
of the operator R z (t, r) is equivalent to a local phase conjugation of the 
signal which enhances the phase derivative related to the Doppler shift 
[30]. Similarly, one usually defines the ambiguity function A z (9,t) as 
the inverse Fourier transform (Ff 1 ) of R z (t, r) with respect to the first 
variable t. 

A z (6,t) = F t - l [z{t + ^)z*(t - T -)\ = F t - l [R z (t,r)} (17) 

Thus, W z (t, v) and A z (9, r) are related by the two-dimensional Fourier 
transform 

W z {t,v) = J J A z {9,T)e-i {te+VT UedT (18) 

The ambiguity function reduces to a deterministic correlation func- 
tion in time if 9 is set to zero. Similarly, as can be seen in a dual 
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form which starts with Z(y) = F t [z(t)\, it can be seen as a de- 
terministic correlation of spectra if r is set to zero. Starting from 
the Wigner-Ville distribution W z (t, u), a wide class of Time- Frequency 
Distributions (TFD) can be developped by means of bi-dimensional 
smoothing operations with appropriate kernels </>(#, r): 



c z (t,v,<t>) = J J J ^- t)6 ~ VT) 4>{e, T)z(i + -) z *{i - -)dtidTde (19) 

The purpose of such a generalisation of the WVD is to provide the 
experimentalist with an operationnal definition of the Time-Frequency 
Distribution TFD z (t, u) of any signal z(t), which can be interpreted as 
a time-frequency power density: TFD z (t,v)AiyAt is the energy of the 
signal z(t) between the instants t and t + At, in the frequency band 
[u v+Au\. Indeed, one can show that with such a definition (smoothing 
kernels), the so-called covariance property (invariance with respect to 
time and frequency translations or modulations) and marginal preser- 
vations are respected : 



PSD z (u) = lim - / TFD z (t,u) dt (20) 

T^oo T J 

-T/2 

+00 

\z(t)z*{t)\ = J TFD z {t,v)du (21) 

—00 

An example of such a Time- Frequency analysis is displayed on figure 
(9) corresponding to the time-serie presented on figure (4). The center 
part of figure (9) is a pseudo-color image (with a linear gray-scale) of 
the level TFD sca t{t,v). This representation evidences the existence of 
finite duration events, with local random Doppler shifts corresponding 
to their own advection velocity. On the left side of figure (9), the time 
averaged spectrum has been represented, together with the frequency 

+T/2 

marginal PSD z {v) = lim ± / TFD scat {t,u) dt of TFD scat {t,u). 

The latter representation demonstrates the blurring effect of the ran- 
dom advection by the large scale velocity field. A way to overcome this 
random sweeping effect is to account for the instantaneous advection 
velocity (at each time where a vorticity event is present, detected by a 
large amplitude of p(t)) . Noting that, as the major effect of the advec- 
tion velocity is to induce a phase modulation according to a Doppler 
frequency shift, we have observed that random advection velocity can 
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Figure 9. Time- Frequency Analysis of z aca t(t) in the LEGI air jet flow (R\ ~ 785). 
Center : TFD z (t,u), Left : PSD marginal (o) and averaged PSD (-), Top : 
instantaneous power marginal 



be safely compensated and fully eliminated by considering only the 
modulus p(t) of the scattering complex pressure signal z scat (t) thus 
ignoring the phase information (j){t). 



4. Time-Space Correlations 

4.1. Characteristic time of a vorticity mode 

Following the later remark, the characteristic time of any spatial vor- 
ticity Fourier mode ^j_(q scat ,f) can be studied by considering only 
its modulus oc p(t). Figure (10) shows the normalized auto-covariance 
Cov pp {t) of the vorticity modulus for different wave- vectors q sca u aligned 
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with the jet axis in the case of the LEGI experiment. The typical 
shape of these correlation curves (right figure) brings out two times: a 
short one (hereafter refered as r s ) characterized by a fast decrease of 
the correlation, and a longer one (noted Tg) associated with low, but 
significant correlation values. Qualitatively, the short time part can be 
rather well fitted by a Gaussian curve (as shown in the left figure), while 
the long time correlation part has rather an exponential-like shape. 




x (ms) x (s) 

Figure 10. Time correlation of vorticity modulus for different scales l/q SC at in the 
LEGI jet (Rx ~ 740) (log-lin coordinates in the insert) 

The short time corresponds to the average time during which the 
component of a vorticity event Q± is tuned with the analyzing wave- 
number q sca tt- This short time t s clearly varies with the analyzing scale 
Qscatt (see below). 

In our experiments, the duration t s was systematically much shorter 
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than the mean transit time in the measurement volume evaluated from 
the size L of the scattering volume and the mean advection velocity: 
t s <C j^— (the mean transit time in this experiment is ~ 25 ms). 

Vavg 

Accordingly, this latter observation reveals the Lagrangian feature of 
the time analysis performed on the spatial Fourier modes of vorticity 
by means of the covariance estimator Cov pp (t). The short time cor- 
responds to the average duration during which the component £l± of 
a vorticity event is observed at the the analyzing wave- number q sca t- 
In our experiments, this duration t s is systematically shorter than the 
corresponding crossing time in the measurement volume. This short 
time t s clearly varies with the analyzing scale q^} a t ■ On the contrary, 
the long time part of the correlation does not seem to depend on the 
scale qjJatf In particular, we observe that the correlation curves cross 
the null value axis at a typical time of the order of tq = 2ir{U rms qo)~ l 
where go is the wave-number characteristic of the integral length scale 
of the flow. In our experimental conditions, this large scale time is 
longer than the crossing time of the measurement volume. This weak 
auto-correlation of the modulus p(t) at long time delays r (larger than 
the duration time r s ) suggests that it should reflect the time behaviour 
of the inter-event dynamic (involving several isolated events). The sig- 
nificantly non zero values, is an indication of the organization of several 
events belonging to the same large structure of the turbulent flow. We 
will not investigate further these long time correlation property as they 
clearly do not correspond to a Lagrangian description of the flow. 
Focusing on the short time r s , the main question is to know how it 
depends on both the length scale and the Reynolds number. 

4.2. Scale dependence 

The figure (11) shows the behaviour of the short time r s with respect 
to the length scale q SC at- Here r s has been arbitrarily defined as the 
width at the half correlation value. We checked that another definition 
of t s (for instance, the standard deviation of the Gaussian fit) does not 
change the results. 

Despite the low extent of the analyzing wave-number range (~ one 
decade), data are in better agreement with a slope —1 than a slope 
—2/3, meaning that the charateristic time of Q±(q S catti t), depends 
linearly on the scale, whatever the Reynolds number. Thus, t s does not 
scale as the eddy turn over time (~ q~ 2 ^ 3 ). For the smallest inertial scale 
range (the only one attainable with our acoustic scattering device), 
experiments show that t s is actually proportional to the local sweeping 
time, defined as T sweep = 2iT(U rms q sca t)~ l ■ Putting forward an analogy 
between a given Fourier mode ^±(q SC att) (in the Fourier space) and a 
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Figure 11. Scaling of r s (normalized by the large scale time to) with the wave 
number q scat : (x) R\ ~ 530), (V) R x ^ 6090. 



velocity increment Sv(5r) over a given separation \deltar ~ q S catt l (in 
the real space) then, the experimental data of the figure (11) are in 
an agreement with those obtained, numerically by Sanada et al. [33], 
and experimentally by Xu et al. [34]. Such a result would suggest that 
Tsweep is the pertinent time of the vortex stretchings and/or tippings 
due to the local velocity gradients. In both papers [33, 34], it appears 
that for the largest scales of the inertial range, the eddy turn-over time 
possibly becomes again shorter than the sweeping one. Unfortunately 
our acoustic device can not reach such inertial scales. 
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4.3. Reynolds dependence 

The figure (11) also indicates that the ratio of the vorticity mode 
life-time r s on the local sweeping time T sweep weakly depends on the 
Reynolds number. 




The figure (12) shows how this ratio a(R\) = — — — behaves with the 
Reynolds number. Thanks to the large range of investigated Reynolds 
numbers, (100 < R\ < 6000), it is reasonable to infer a power law 
scaling with an exponent close to —1/4. Note that the numerical data 
of Sanada et al. [33] are in very good agreement with our experiments 
(the two points of figure 12 have been taken from the figure 8 of [33]). 
In others words, for a fixed scale 2-7r(q sca t)~ 1 , the vorticity correla- 
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tion time decreases with the Reynolds number. In connection with the 
vorticity time-correlation, Novikov [35, 36] has introduced a typical 
characteristic length £ c which is the inertial scale where the viscous 
dissipation of the enstrophy is balanced by the large scale effects. In 

—2/5 

terms of time scales, the Novikov prediction is such as r c ~ tq1A0R x 
(where r c is the time scale of the vorticity correlation). 
This prediction can not be directly compared to the experimental scal- 
ing of the figure (12) because the latter relies on the analysis of a single 
length scale whereas the theoretical result is based on the global bud- 
get of the vorticity correlations involving all the scales. However, one 
expects that experimental data are not in agreement with the previous 
prediction, because it is mainly derived from Kolmogorov arguments 
and it does not account for sweeping effects. 



5. Lagrangian Velocity Measurements 

The interest in the statistics of the Lagrangian velocity (following the 
fluid particules) in turbulent flows has been recently renewed by the 
emergence of new experimental techniques for the continuous tracking 
along time of isolated particules [37, 38] (using acoustical devices) or 
the measurement of a large number of individual particule acceleration 
events [39] (using fast optical techniques). In the first type of experi- 
ment [37], a very few number (ideally one) of solid beads are injected 
in a closed turbulent water flow. The trajectory of a single particule is 
then followed along time intervals, as large as the integral length scale 
of the turbulent flow, using an acoustic scattering device consisting in 
an array of small transducers (in a backscattering configuration). 
The purpose of the present experiment is to extend the latter acoustic 
tracking technique to the study of the Lagrangian turbulent velocity 
in open turbulent air flows. Material particules, with an acoustic index 
different from the acoustic index of the propagating medium, can effi- 
ciently scatter acoustic waves [4]. Provided the acoustic contrast (the 
relative difference between the acoustic index of the particule and that 
of the medium) is sufficient, acoustic scattering by particules injected 
in a turbulent flow can largely overcome the acoustic scattering by the 
turbulent vorticity fluctuations. When a particule is advected by the 
local flow velocity, the scattered pressure amplitude p sca t(t), at a large 
enough scattering wave-number {q sca t r > 2-7T, where r is the typical 
size of the particule), will exhibit Doppler shift phase modulations. We 
have seeded the LEGI jet, with very small soap bubbles (of diameter 
< 3 mm), at a small rate of injection (typ. < 100 bubbles/ s) to have at 
most one particule at the same time in the scaterring volume V sca t- 
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The soap bubbles are inflated with gazeous helium, with a home-made 
apparatus, so as to match pecisely the average density of the bubble 
with the density of the surrounding fluid (the low density helium com- 
pensates properly the mass of the soap film). Thanks to this density 
matching, and to the small size of the bubbles (limiting wake effects 
behind the bubble [37]), inertial effects are small and each isolated 
particule is expected to follow, nearly instantaneously, the fluid par- 
ticule velocity. In this experiment, the scattering angle is set to 120° 
in order to increase, both the wave-number q S cat 

and the size of the 

scatterring volume V sca t in the direction of q sca t (see figure (3)) up 
to a length larger than the integral length scale Li nt of the turbulent 
flow. Moreover, at such a large scattering angle, the contribution of 
the vorticity fluctuations to the scattering amplitude is significantly 
reduced (see figure (2)). 

In this preliminary experiment, the mean velocity was about 4 m/s 
(R x ~ 530). 

5.1. Scattered pressure signal PSD and covariance 

A first experiment has been performed to study the correlation and 
the spectrum of the acoustic pressure scattered by the soap bubbles. 
Simultaneous acquisitions have been realized with the superposition of 
two incident frequencies (80 kHz — 100 kHz and 100 kHz — 120 kHz). 
The two demodulated scattered signals around each incident frequency 
is then extracted by means of a band pass filtering and a demodulation. 
Such a double channel experiment allows a direct comparison of the de- 
tection at different incoming sound frequencies. In all experiments, the 
passage of a bubble in the scattering volume is simultaneously detected 
(around each v Q ) as a strong increase of the scattered pressure intensity 
I SC at(t) = p{t) 2 -> indicating a large signal to noise ratio of the acoustic 
scattering set-up. Also, the Doppler shift frequency, around each incom- 
ing frequency v scales continuously in time with v Q according to Eqn. 
12. The high level of similarity at each time t, of both p(t) and ^jp- is 
a good indication of the efficiency of the acoustic scattering technique. 
Moreover, thanks to the small bubble flow rate, we effectively observe 
isolated bubbles passing through the scattering volume. 
Indeed, the time interval during which the scattered pressure is signif- 
icant is of order AT = x^-, where AL is the size of the scattering 

Vavg 

volume in the direction of q sca t (aligned with the mean flow). This 
is illustrated by the plot of the normalized auto-covariance of p(t) on 
figure (13-left). The covariance, decreases slowly from its maximum 
value (at the delay r = 0), to a null value at the delay r = ~ 50 ms. 
Remind that in the case of acoustic scattering by vorticity fluctuations, 
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Figure 13. Normalized Covariance C pp (t) of p sca t(i) (left) and PSD aca t(5v) (right) 
of the acoustic pressure scattered by soap bubbles in the LEGI turbulent jet flow at 
Rx ~ 500. 

a much more rapid decrease of the auto-covariance is observed leading 
to an estimation of the duration of vorticity events of a few ms (cf 
figure (10)). The length of the negative part of the auto-covariance can 
be ascribed to the mean time between the passage of two consecutive 
bubbles in the scattering volume. The demodulated scattered pressure 
signals have been largely over-sampled at F s = 32768 Hz and further 
band-pass filtered with a Butterworth digital filter of order 8. The right 
part of figure (13) displays the power spectral density for the scattered 
pressure signal at an incoming frequency v = 80 kHz. As expected, 
the shape of PSD scat {5v) is Gaussian (8v avg and 5v st d, in excellent 
agreement with the statistics of the longitudinal flow velocity (given by 
hot-wire anemometry). 
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5.2. Lagrangian velocity PDF AND spectrum 

To investigate the statistical properties of the Lagrangian velocity, a 
very large number of bubble trajectories is required for purpose of 
statistical convergency. 
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Figure U. PDF (left) and PSD (right) of the Lagrangian Velocity in the LEGI 
turbulent jet flow at R\ ~ 500 

Thus we have performed a long time experiment, collecting 2.7 10 s 
complex data points sampled at F s = 32768 Hz. Using a post- 
processing algorithm, relying on an appropriate thresholding of the 
scattered acoustic intensity I sca t(t), the passage of isolated bubbles in 
the scattering volume are identified and localized. 
About 13 10 6 passing bubbles have been identified and localized. After 
a band-pass filtering, the projection of the instantaneous velocity (onto 
the direction of the jet axis) of each isolated bubble is estimated from 
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the unwrapped phase signal <j)(t) (according to Vb u bble(t) = ^— t ^^-) 
with a finite difference schema. 

From all these Lagrangian velocity signal we have computed firstly the 
probability density function PDF(vb u bbie) of the Lagrangian velocity 
depicted on figure (14 left). The Lagrangian PDF has a nice gaussian 
shape, indicating a mean value and rms values of the Lagrangian veloc- 
ity (note that they are ensemble averages) close to the corresponding 
Eulerian values (determined with hot-wire anemometry). The equality 
of the Lagrangian and Eulerian probability densities is a well known 
property of homogeneous turbulent flows [2, 40]. 

The time power density spectrum PSDi, a g(5v) have been computed 
with the Welch's averaged periodogram method. The resulting PSD 
is displayed on figure (14 right). In this log-log representation, a clear 
power-law scaling over roughly one decade is visible. The exponent of 
the observed power law scaling is closed to the value —2 in agreement 
with the theoretical prediction of Kolmogorov [2]. At higher frequencies, 
the DSP display a significant cut-off, above the frequency v c ~ 70 Hz 
very close to the low-pass filtering frequency due to the radius rf, = 
3 mm of the bubbles (z/& ~ e 1 / 3 r b 2 ^ 3 ~ 70 Hz where e is the mean 
energy dissipation rate). 

6. Concluding remarks 

Wave scattering (e.g. light and neutron scattering in condensed matter 
and X-ray scattering in solid state physics) has proved to be a powerful 
tool in the study and of complex systems (e.g. phase transitions, DNA 
structure and denaturation) . Turbulent flows are also known to be very 
complex and disordered. The scattering of coherent acoustic waves by 
turbulent vorticity fluctuations allows the direct and continuous (in 
time) probing of a spatial Fourier mode Q±(q scat , t), at a well defined 
spatial wave- vector q sca f. In addition, the acoustic scattering technique 
is non intrusive and sensitive to the local orientation of the vorticity 
field [41]. Thus it appears as a relevant means of investigation of turbu- 
lence, where the vorticity dynamic is known to play a crucial role (e.g. 
vortex stretching and bending), particularly in the energy transfers 
across the length scales [43, 44]. By selected different spatial wave- 
vectors, one is able to analyse scale by scale the statistical properties 
of the turbulent vorticity distribution. For example, we have evidenced 
a scaling law with exponent 1/3 for enstrophy spectral density, in 
agreement with the Kolmogorov prediction. Time correlation of the 
amplitude (modulus) vorticity modes reveals a short life-time of the 
vorticity events at a fixed scale (possibly coherent structures [32]). As 
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the average duration of the detected vorticity events is much smaller 
than their mean transit time across the scattering volume we argue that 
our analysis is of a Lagrangian type. Still, the vorticity correlations 
also exhibit a long time behaviour, with a lower level of correlation, 
indicating probably a large scale structure of the vorticity distribution 
which could be interpreted as the consequence of the turbulent cascade 
process. 

Acoustic scattering can also be used to access Lagrangian statistics of 
the turbulent velocity thanks to the possibility to detect very small 
particules (using high enough ultrasound frequencies), on very large 
volumes (of order the integral length scale of the flow) . With this pre- 
liminary experiment, we have found again, in an open flow, Lagrangian 
probability density function of the turbulent velocity field equal to the 
Eulerian PDF as expected. Beside, we have also evidenced a Lagrangian 
power law scaling (oc lo~ 2 ) for spectral power density, according to the 
Kolmogorov 41 theory. The power law scaling extends over roughly 
one decade. Further investigation are in progress. Inparticular, we have 
already generated and successfully detected bubbles with a smaller 
radius (< 1 mm) and thus we expect to increase the latter power law 
scaling range towards the dissipative time-scales. 
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